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ABSTRACT: Thermally conductive resins are needed for 
bipolar plates in fuel cells. Currently, the materials used for 
these bipolar plates often contain a single type of graphite in 
a thermosetting resin, m this study, varying amounts of four 
different types of polyacrylonitrile carbon fillers (Ketjenblack 
carbon black, Thermocarb synthetic graphite, Fortafil 243 car- 
bon fiber, and Panex 30 carbon fiber) were added to a ther- 
moplastic matrix (Vectra A950RX Liquid Crystal Polymer), 
with die resulting resins tested for through-plane and in- 
plane thermal conductivity. There are two unique contribu- 
tions of this work. The first contribution is the use of the 
Nielsen model for die through-plane thermal conductivity as 



a function of the single filler volume fraction. The model fits 
the data for all composites well. The second contribution is 
the development of a new, accurate, empirical model to pre- 
dict the in-plane thermal conductivity for all resins contain- 
ing synthetic graphite or carbon fiber. Both of these models 
will form the basis for the development of new thermal con- 
ductivity models for composites with multiple fillers for fuel 
cell bipolar plate applications. © 2007 Wiley Periodicals, Inc. 
J Appl Polym Sd 105: 3309-3316, 2007 
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INTRODUCTION 

Most polymer resins are thermally insulating. 
Increasing the thermal conductivity of these resins 
allows them to be used in other applications. One 
emerging market for thermally conductive resins is 
for bipolar plates for use in fuel cells. The bipolar 
plate separates one cell from the next, with the plate 
carrying hydrogen gas on one side and air (oxygen) 
on the other side: Bipolar plates require high thermal 
conductivity (to conduct heat), low gas permeability, 
and good dimensional stability. 

Typical thermal conductivity values for some com- 
mon materials are 0.2-0.3 for polymers, 234 for alumi- 
num, 400 for copper, and 600 for graphite (all values 
in W m -1 K -1 ). One approach to improve the ther- 
mal conductivity of a polymer is through the addition 
of a conductive filler material such as carbon and 
metal. Conductive resins with a thermal conductivity 
from ~1 to 30 W m" 1 KT 1 can be used in heat sink 
applications. 1 
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A significant amount of work has been conducted 
varying the amount of single conductive fillers in a 
composite material. 2 " 9 For example, ceramic fibers/ 
particles (boron nitride, aluminum nitride, alumi- 
num oxide), metal fibers/particles (aluminum, steel, 
iron, copper, silver), and Ni-coated glass fibers have 
been used. 3 ' 10 " 13 Metallic fillers have several disad- 
vantages, relative to carbon, which include higher 
density and greater susceptibility to oxidation. Vari- 
ous types , of .carbons have been effective^ conductive 
fi]|erd._ Pgr . e^uiipler'Syndietic graphite particles and 
carbon fibers are often added to jpolyinere to increase 
the composite thermal conductivity. W 4 ' 15 

In this project, researchers performed compound- 
ing runs followed by injection molding of carbon 
fibef/Hquid ^crystal "polymer (LCP) ."test ~ specimens. 
Varying amounts of four different types of carbon 
fillers (Ketjenblack carbon black, Thermocarb syn- 
thetic graphite, Fortafil 243 carbon fiber, and Panex 
30 carbon. fiber) were; added to Vectra A950RX LCP.. 
The resulting. "Single "filler composites were then 
tested for thermal conductivity. 

We note that a long term goal of this research pro- - 
ject is to understand the effects of multiple carbon 
fillers on the mermal properties of fuel cell bipqlar 
plates . Fundamental research is required to attain 
this goal, namely me experimental measurement and 
model formulation, for the single filler composites, 
with regard to ih-plahe and through-plane thermal 
conductivities. These results are presented here. For 
• experiments use the, guardecl heat;flow 
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TABLE I 

Properties of Ticona's Vectra A950RX 34 



Melting point 
Tensile modulus 

(1 irun/mm) 
Tensile stress at break 

(5 mm/min) 
Tensile strain at break 

(5 mm/min) 
Flexural modulus at 23°C 
Notched izod impact 

strength at 23°C 
Density at 23°C 
Volumetric electrical 

conductivity at 23°C 
Surface electrical conductivity 
Thermal conductivity at 23°C 
Humidity absorption 

(23°C/50% RH) 
Mold shrinkage-parallel 
Mold shrinkage-normal 
Coefficient of linear thermal 

expansion-parallel 
Coefficient of linear thermal 

exp ansion-normal 



280°C 

10.6 GPa 

182 MPa 

3.4% 
9.1 GPa 

95 kj/m 2 
1.40 g/cm 3 

10" 15 S/cm 
10 -14 S 
0.2 W m" 1 K" 1 (approx.) 

0.03 wt % 
0.0% 
0.7% 

0.04 X 10 - V°C 

0.38 X 10 -4 /°C 



method to measure through-plane thermal conduc- 
tivity 16 and the transient plane source method/ 7 " 21 a 
relatively new analytical technique, which allows for 
measurement of both the in-plane and through-plane 
thermal conductivities in a single experiment. 

There are two specific modeling-related goals of 
this research, which are new contributions to the 
composites literature. The first goal was to apply 
Nielsen's thermal conductivity model 22 for the 
through-plane thermal conductivity of these formu- 
lations. The second goal is to use a simple empirical 
model developed earlier in our group, ' along 
with the Nielsen model results, to predict ih-plane 
thermal conductivities. We note that the technical 
literature has extensive experimental and modeling 
data for through-plane thermal conductivity, 22 - 25 " 
but there is little reported experimental data nor 
model development for the in-plane thermal con- 
ductivity. 

MATERIALS AND EXPERIMENTAL METHODS 
Materials 

The matrix used for this project was Ticona's 
(Summit, NJ) Vectra A950RX LCP, which is a highly 
ordered thermoplastic copolymer consisting of 73 
mol % hydroxybenzoic acid and 27 mol % hydroxy- 
naphtholic acid, this LCP has the properties needed 
for bipolar plates, namely high-dimensional stability 
up to a . temperature of 250° C, extremely short mold- 
ing times (often 5-10 s), exceptional dimensional 
reproducibility, ^ chemically resistant in acidic envi- 
ronments present in a fuel cell, and a low hydrogen 



gas permeation rate. 3435 The properties of this poly- 
mer are shown in Table I. 34 

The first carbon filler used in this study was Ket- 
jeriblack EC-600 JD, carbon black This is an electri- 
cally conductive carbon black available from Akzo 
Nobel (Chicago, IL). The highly branched, high sur- 
face area carbon black structure allows it to contact a 
large amount of polymer, which results in improved 
electrical conductivity at low carbon black concentra- 
tions (often 5-7 wt %). The properties of carbon 
black are given in Table n. 36 According to the ven- 
dor literature, carbon black is. sold in the form of 
pellets that are 100 um-2 mm in size and, upon mix- 
ing into a polymer, easily separates into primary 
aggregates 30-100-nm long. 36 A diagram of the car- 
bon black structure is shown in vendor literature 36 
and prior work from our group 37 It is noted that the 
thermal conductivity of pure carbon black cannot be 
estimated because of the small size of the aggre- 
gates. 

Table HI shows the properties of the Asbury Car- 
bons' (Asbury, NJ) Thermocarb TC-300, which is a 
primary synthetic graphite that was previously sold 
by Conoco (Houston, TX). 38 ' 39 Thermocarb TC-300 
has high electrical and thermal conductivity, and a 
mean particle aspect ratio of ~ 1.7. Thermocarb TC- 
300 is produced from a thermally treated, highly aro- 
matic petroleum feedstock, and contains very few 
impurities. A photomicrograph of this synthetic 
graphite is shown in the vendor literature 38 and 
prior work from our group 37 

Fortafil 243, sold by Toho Tenax America (Rock- 
wood, NJ) is a pplyacrylonitrile (PAN)-based 3.2-mm 
chopped and pelletized carbon fiber used to improve 
the electrical and thermal conductivity of the resin. 
Fortafil 243 was surface-treated and then formed 
into pellets. A proprietary polymer (sizing) is used 
as a binder for the pellets that also promotes adhe- 
sion with the matrix. Table IV shows the properties 
of this carbon fiber, which is 95 wt % carbon. 

Table V shows the properties of Zoltek's (St. 
Louis, MO) Panex 30 MF milled 150-um long high 
purity carbon fiber. This carbon fiber is PAN-based 
and is electrochemically surface-treated but not 
sized. Panex 30 is produced by a high-temperature 



TABLE n 

Properties of Akzo Nobel Ketjenblack EC 600-JD 

Carbon Black 36 



Electrical resistivity 
Aggregate size 
Specific gravity 
Apparent bulk density 
Ash content, max 
Moisture, max 

Brunauer-Emmett-Teller surface area 
Pore volume 



0.01-0.1 Q cm 
30-100 nm 
1.8 g/cm 3 
100-120 kg/m 3 
0.1 wt % 
0.5 wt % 
1250 m 2 /g 
480-510 cm 3 /100 g 
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TABLE m 



Properties of Thermocarb TC-300 Synthetic Graphite 



38,39 



Filler 



Thermocarb TC-300 
synthetic graphite 



Carbon content (wt %) 

Ash (wt %) 

Sulfur (wt %) 

Density (g/cm 3 ) 

B runauer-Emme tt-Tell er 

surface area (m 2 /g) 
Thermal conductivity 

at 23°C (W irT 1 K" 1 ) 

Electrical conductivity 

of bulk carbon powder 

at 150 i|f, 23°C, parallel 

to pressing axis (S/cm) 
Particle shape 
Particle aspect ratio 
Sieve analysis, microns (wt %) 

+600 um 

+500 um 

+300 um 

+212 um 

+ 180 um 

+ 150 um 

+75 um 

+44 um 

-44 um 



99.91 
<0.1 

0.004 

2.24 

1.4 

600 in. "a" 
crystallographic 
direction 



50 
Acicular 
1.7 

0.19 
0.36 
5.24 
12.04 
8.25 
12.44 
34.89 
16.17 
10.42 



batch graphitization process that produces fiber that 
is 99.5 wt % carbon. Panex 30 also has a high elec- 
trical and thermal conductivity. 

Thermal conductivity was measured on compo- 
sites containing varying amounts of these carbon 
fibers in Vectra A950RX. The concentrations. (shown 
in wt % and the corresponding vol %) for these sin- 
gle filler composites are shown in Table VI. 

Test specimen fabrication 

For this entire project, the fillers were used as- 
received. Vectra A950RX was dried in an indirect- 
heated dehumidifying drying oven at 150°C and 
then stored in moisture barrier bags. 



TABLE IV 



Properties of Fortafil 243 Carbon Fiber 



.40 



Carbon content 
Electrical resistivity 
Thermal conductivity 
Tensile strength 
Tensile modulus 
Specific gravity 
Fiber diameter 
Fiber shape 
Fiber mean length 
Binder content 



Bulk density 



95 wt % 
0.00167 0 cm 

20 W m" 1 K" 1 (axial direction) 

3800 MPa 

227 GPa 

1.74 g/cm 3 

7.3 um 

Round 

3.2 mm (entire range is 2.3-4.1 mm) 
2.6 wt % proprietary polymer 

that adheres pellet together and 

promotes adhesion with 

nylon matrix 
356 g/L 



TABLE V 

Properties of Panex 30 MF Milled High Purity 

Carbon Fiber 41 



Carbon content 
Electrical resistivity 
Thermal conductivity 

Tensile strength 
Tensile modulus 
Specific gravity 
Fiber diameter 
Fiber shape 
Fiber mean length 
Bulk density 



99.5 wt % 
0.0014 fi cm 
22 W rrT 1 K" 1 (axial. direction, 
approximate) 
3600 MPa 
207 GPa 
1.75 g/cm 3 
7.4 um 
Round 
150 um 
445 g/L 



The extruder used was an American Leistritz Ex- 
truder Corp. (Somerville, NJ) model ZSE 27. This ex- 
truder has a 27-mm corotating mtermeshing twin 
screw with 10 zones and a length/diameter ratio of 
40. The screw design used is the same as that used 
in prior published research. 42 The screw design was 
chosen to allow a large concentration of filler to mix 
with the matrix material and thereby achieve the 
maximum possible conductivity. The Vectra polymer 
pellets were introduced in Zone 1. A side stuff er 
located at Zone 5 was used to introduce the carbon 
fillers into the polymer melt. Two Schenck (White- 
water, Wl) AccuRate gravimetric feeders were used 
to accurately control the amount of each material 
added to the extruder. 

After passing through the extruder, the composite 
strands (3 mm in diameter) entered a water bath 
and then a pelletizer that produced nominally 3-mm 
long pellets. After compounding, the pelletized com- 



TABLE VI 
Single Filler Loading Levels 



in Vectra A950RX 





Carbon 


Thermocarb 






Filler 


black 


TC-300 


Fortafil 243 


Panex 30 


(wt %) 


(vol %) 


(vol %) 


(vol%) 


(vol %) 


2.5 


1.9 


N/A 


N/A 


N/A 


4.0 


3.1 


N/A 


N/A 


N/A 


5.0 


3.9 


N/A 


4.1 


4.0 


6.0 


4.7 


N/A 


N/A 


N/A 


7.5 


6.0 


N/A 


6.1 


6.1 


10 


8.0 


6.5 


8.2 


8.2 


15 


12.1 


9.9 


12.4 


12.4 


20 


N/A 


13.5 


16.8 


16.7 


25 


N/A 


17.2 


21.2 


21.1 


30 


• N/A . 


21.1 


25.5 


25.5 


35 


N/A 


25.2 


30.2 


30.1 


40 


N/A 


29.3 


34.9 


34.8 


45 


-N/A 


33.8 


39.7 


39.9 


50 


N/A 


38.5 


44.6 


44.4 


55 


N/A 


43.3 


49.6 


49.4 


60 


N/A 


48.4 


54.7 


54.5 


65 


' N/A 


53.7 


N/A 


N/A 


70 


N/A 


59.3 


N/A 


N/A 


75 


N/A 


65.2 


N/A 


N/A 
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posite resin was dried again and then stored in 
moisture barrier bags prior to injection molding. 

A Niigata (Itasca, IL) injection molding machine, 
model NE85UA4, was used to produce test speci- 
mens. This machine has a 40-mm diameter single 
screw with a length/diameter ratio of 18. The lengths 
of the feed, compression, and metering sections of 
the single screw are 396, 180, and 144 mm, respec- 
tively. A four-cavity mold was used to produce 6.4- 
cm diameter (end gated) 3.2-mm thick disks. The 
thermal conductivity of all formulations was deter- 
mined. Prior to conducting thermal conductivity 
tests, the samples were conditioned at 23° C and 50% 
relative humidity for 88 h and then tested. 43 

Filler length and orientation test method 

Because of the small size of the carbon black (the pri- 
mary aggregates were 30-100 nm), the carbon black 
filler length and orientation were not measured. The 
following procedure was used to determine the 
length of the other carbon fillers in the molded test 
specimens. 33 ' 44 Diemylmetriamine was used to dis- 
solve the matrix. The fillers were then dispersed onto 
a glass slide and viewed using an Olympus (Center 
Valley, PA) SZH10 reflected light microscope at a 
magnification of 60X. Digital images were taken of 
the filler and were then analyzed with Adobe (San 
Jose, CA) Photoshop version 5.0 installed with the 
Image Processing Toolkit. 

The following procedure was used to determine 
the orientation of the synthetic graphite and carbon 
fibers 33/44 For the thermal conductivity samples, the 
center portion was cut out of a disk and set in epoxy 
such that the through sample thickness (3.2 mm) 
face could be viewed. The samples were then pol- 
ished and viewed using an Olympus BX60 reflected 
light microscope at a magnification of 200 X. 

Through-plane thermal conductivity test method 

The through-plane thermal conductivity of a 3.2-mm 
thick, 5-cm diameter disc-shaped test specimen was 
measured at 55°C using a Holometrix (Burlington, 
MA) Model TCA-300 Thermal Conductivity Ana- 
lyzer according to the ASTM F433 guarded heat flow 
meter method. 16 For each formulation, six samples 
were tested. 

Transient plane source thermal conductivity 
test method 

The Mathis Instruments (Halifax, Nova Scotia, Can- 
ada) Hot Disk Thermal Constants Analyzer is an 
emerging technology that can measure the in-plane 
and through-plane thermal conductivity of an aniso- 
tropic material in the same test, using the transient 
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plane source technique. 17 " 24 ' 45 ' 46 The sensor used in 
this test method consisted of a 10-um thick nickel 
foil embedded between two 25.4-um-thick layers of 
Kapton polyimide film. The nickel foil was wound 
in a double spiral pattern with a radius R of 3.189 
mm. This sensor is placed between two polymer com- 
posite disks of diameter D = 63.5 mm and thickness T 
= 3.18 mm. The thermal conductivities were meas- 
ured at 23°C. During the experiment, the sensor was 
heated by a constant electrical current (about 0.05 
W) over a short period of time (typically 5 s). The 
generated heat dissipated within the double spiral is 
conducted through the kapton insulating layer and 
into the surrounding sample, causing a rise in the 
temperature of the sensor and the sample. The average 
transient temperature increase of the sensor is simulta- 
neously measured by recording the change in electri- 
cal resistance of the nickel sensor and can then be used 
to back-calculate the in-plane and through-plane ther- 
mal conductivities of the polymer composite. 

RESULTS 

Filler length and orientation results 

For the molded test specimens containing synthetic 
graphite, the filler length and aspect ratio were typi- 
cally 50 um and 1.68, respectively. The synthetic 
graphite length results are comparable to that of 
Heiser and King 33 For the molded test specimens 
containing both carbon fibers, the fiber length was 
typically 70 um. The corresponding fiber aspect ratio 
(length/diameter) was 9. The carbon fiber length 
results are comparable to that of Heiser and King 33 
and Bigg 47 

As shown previously for me through-plane ther- 
mal conductivity samples, the fillers are primarily 
oriented transverse to the conductivity measurement 
direction. These results are comparable to our prior 
work on filler orientation 33 Photomicrographs are 
available in the published literature. 33 ' 42 

Thermal conductivity experimental results 

The filler size, shape, concentration, dispersion (degree 
of mixing), orientation, bonding between, the filler and 
matrix, thermal conductivity of the constituents (filler 
and matrix), and the crystallinity of the polymer 
(increasing crystallinity improves merrrial conductiv- 
ity) are factors that affect the thermal conductivity of a 
composite. Figures 1 and 2 display the experimental 
through-plane thermal conductivity data ^through * or 
the samples containing Ketjenblack, Therrriocarb, For- 
tafil 243, and Panex 30. Several generalizations can be 
made regarding these data. First of all, it is noted that 
at higher loading levels, the Thefmocarb composites ex- 
hibited the highest mermal ccniductryity of the filled 
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0 5r 




0.4 




0.3 


1 






0.2^ 




0.1 



°0 0.025 0.05 0.075 0.1 0.125 0.15 

4, Volume Fraction 

Figure 1 Through-plane experimental (data points)^ and 
theoretical (line) thermal conductivities in W m" 1 K 1 for 
composites containing Ketjenblack (diamonds). 

composites tested. This can be attributed to the high ther- 
mal conductivity of the Thermocarb (600 W m" 1 K" 1 ) 
when compared with the carbon fibers (both about 
20 W m _1 K" 1 ) and carbon black (thermal conductiv- 
ity measurement was not available, but we estimate 
it as 2.1 W m" 1 K" 1 in the modeling section of this 
article). Second, the Ketjenblack composites exhibit a 
relatively linear increase in through-plane thermal 
conductivity over the range of filler volume fractions 
tested in this research. We could not extrude or injec- 
tion mold composites with filler loadings above 12.1 
vol % due to the large increase in composite viscosity, 
as we note in a prior publication. 42 We note that for 
the synthetic graphite or carbon fiber composites, the 
thermal conductivity is relatively linear for volume 
fractions less than 0.12, after which there is a rapid 




°0 0.1 0.2 0.3 0.4 0.5 0.6 0:7 
4>, Volume Fraction 



Figure 2 Through-plane experimental (data points) and 
theoretical (lines) thermal conductivities in W m -1 K _1 for 
composites containing Thermocarb (squares), Fortafil 243 
(triangles), and Panex 30 (circles). - : 




<|>, Volume Fraction 



Figure 3 In-plane experimental (data points) and theoreti- 
cal (line) thermal conductivities in W m" 1 K" 1 for compo- 
sites containing Thermocarb (squares). 

growth. This rapid growth can likely be attributed to 
significant conductive pathways formed in the com- 
posite. 

The in-plane thermal conductivities, for Thermo- 
carb, Fortafil 243, and Panex 30 are shown as a function 
of filler volume fraction $ in Figures 3 and 4. Compo- 
sites containing Ketjenblack were not tested for in-plane 
thermal conductivity, as their thermal conductivity is 
too low for the test method considered here. Composites 
ccmtairung Thermocarb have higher in-plane thermal 
ccmducuvities for all filler loading levels. It is noted that 
the Thermocarb/Vectra composites likely have higher 
in-plane thermal conductivity because of the higher 
thermal conductivity of the constituent Thermocarb (600 
W m" 1 K" 1 ) compared with the Fortafil 243 and Panex 
30 fibers (both about 20 W m" 1 K" 1 ). 

5 1 1 1 — » « ' 



4.5 




°0 0.1 0.2 0.3 0.4 0.5 0.6 
<(>, Volume Fraction 

Figure 4 In-plane experimental (data points) and theoreti- 
cal (line) thermal conductivities in W m _1 K" 1 for compo- 
sites containing Fortafil 243 (triangles) and Panex 30 
(circles). ' 
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Thermal conductivity modeling background 

In many engineering applications, including that for 
bipolar plates' in fuel cells; composite materials are 
thin with respect to one dimension. Many of the ex- 
perimental techniques, and therefore development of 
mermal conductivity models, have been applied to- 
ward predicting the through-plane thermal conduc- 
tivity of a composite. However, as new experimental 
techniques have been developed to measure the in- 
plane thermal conductivity, 21 it is also desired to 
develop models for this transport property as well. 
After briefly reviewing the modeling literature, we 
describe our modeling approach for the anisotropic 
thermal conductivities of the composites tested here. 

For simplicity as well as direct application to the 
research presented here, we will focus on a system 
with one filler. The following nomenclature will be 
used: K represents the thermal conductivity of the 
composite, fcf the thermal conductivity of the ith con- 
stituent (1 = matrix, 2 = filler), and <|) the volume 
fraction of the filler (note 1 - <|> is the matrix volume 
fraction). 

The most basic thermal-conductivity models are 
the standard rule of mixtures [eq. (1)], inverse rule 
of mixtures [eq. (2)], and the geometric rule of mix- 
tures [eq. (3)]. 28 



K = (1 - <t>)fci + #2 

i _ (i - 40 . 4> 



(1) 



(2) 



(3) 



It is noted that these models work well for continu- 
ous, periodic geometries, such as stratified slabs in a 
parallel (rule of mixtures) or transverse (inverse rule 
of mixtures) orientation with the direction of heat 
conduction. However, if the materials are not contin- 
uous, as is the case for carbon-filled resins, the mod- 
els are not accurate. 

There are many advanced models for composite 
systems. The basis of many of them is from Maxwell 
and is based upon potential theory to obtain an exact 
solution for the conductivity of a system with spheri- 
cal noninteracting particles in a continuous matrix. 25 
As most systems do not have noninteracting spheri- 
cal particles, this model is restricted in its utility. 
Hamilton and Crosser 27 developed a semitheoretical 
model that can account for filler shape. The parame- 
ter \|/ indicates the sphericity of the particles, typi- 
cally between 0.58 and 1.0 for the investigated data 
sets. Their model is given as 

r . / >2+(H-l)ti + (n-l)»fe-ti)\ r4rf 
K \ * 2 +(n-l)*i-<K*2-*i) ) 



3 

n = T 



(4b) 



Wh en n = 3, their mqdelxeduces to the Maxwell 
model. The Bruggeman 26 theoretical model also fol- 
lows Maxwell's derivation, but with different 
assumptions about the permeability and field 
strength of the system. It is given as 



1/3 



(5) 



Other advanced models can be derived from the hy- 
pothesis that the ratio of thermal conductivity of a 
composite material to that of the pure polymer is 
proportional to the increase in viscosity, tensile mod- 
ulus, and /or shear modulus. The original basis for 
many of these models is Albert Einstein's model for 
the viscosity of a fluid with dispersed spheres. On 
the basis of this premise, McCullough and cow- 
orkers 29 ' 30 proposed a generalized equation in com- 
bination with traditional mixing rules and a refer- 
ence state, used in composites that show orthotropic 
symmetry. The McCullough model, with some rear- 
rangement, serves as the basis for the Halpin-Tsai 
equations 28 : 



K = 1 + 4X* 
k\ 1 — X<f> 



(6a) 



X = 



k %zl 



(6b) 



where £ varies depending upon the filler used. 
Lewis and Nielsen 22 improved upon the Halpin-Tsai 
model by adding a value to the denominator of the 
primary equation, given by the symbol ^ to take 
into account the filler shape, orientation, and the 
packing. The model is given as, 



K 1 -hAB<t> 



(7a) 



B = 



\(/ ^ 1 + 



— 5— 4> 



(7b) 



(7c) 



m 



where <t> m term is the maximum volumetric packing 
fraction of the filler. The parameters A and <j> m has 
been tabulated for some filler types and orientation 
in the literature. 13 
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Progelhof et al. reviewed many composite ther- 
mal-conductivity models. This review determined 
that the Nielsen model fits the data the best over the 
given data range (0-30 vol %) for a two-phase sys- 
tem. This review studied 62-88-um diameter glass 
spheres and 62-125-um diameter magnesium oxide 
powder in polyethylene. Work done by Weber and 
colleagues 31 ' 32 has expanded the Nielsen model into 
a form more applicable for multiple-filler systems. 
The primary change in the system of equations 
occurs in eq. (7a), where the contributions from fill- 
ers are accounted for in a summation term. 

Recently, with the advent of the transient plane 
source method, our group has been able to obtain 
in-plane and through-plane thermal conductivity 
data. An empirical correlation 23 ' 24 has been devel- 
oped for the square root of the product of the in- 
plane and through-plane thermal conductivities 
given as 

^ fenfehrough = C e D * (8) 
where C and D are adjustable parameters. 

Through-plane thermal conductivity 
modelmg results 

In this section we use Nielsen's model to predict the 
through-plane thermal conductivity fehiough. of con- 
ductive resins in (W m" 1 K -1 ) as a function of the 
filler volume fraction (J). It is noted that Nielsen's 
model typically uses fixed data for the parameters A 
and <$> m , depending on the filler shape, aspect ratio, 
and packing. However, these models have been 
shown either to underestimate thermal conductivity 
or break down at higher filler concentra- 
tions 2 ' 31 ' 32 ' 46 ' 48 Thus, the theoretical through-plane 
thermal conductivities fehrough for Thermocarb, Forta- 
fil 243, and Panex were calculated using eq. (7) by 
adjusting the values of the shape factor A and maxi- 
mum packing fraction <\> m . In all cases, fci — 0.2 W 
m -i K -i (j a bi e rj) The results are shown in Figure 

2 with parameters given as follows: 

Thermocarb: k 2 = 600 W m" 1 K" 1 (Table HI), A = 

6.0, (\> m = 1.00 
Fortafil 243: k 2 = 20 W m" 1 K (Table IV), A = 1.7, 

<l>m = 0.80 

Panex 30: k 2 = 22 W m" 1 KT 1 (Table V), A = 2.0, 
<t>m = 0.72 



bon black is variable, and was used here as an ad- 
justable parameter along with the shape factor. The 
results of the model are shown in Figure 1 with the 
parameters given as 

Carbon black: k 2 = 2.1 W m" 1 FT 1 , A = 1720, (j> m 
= 0.64 

In-plane thermal conductivity modeling results 

It was also desired as part of this work to develop a 
simple model for the in-plane thermal conductivity 
fen by combining eqs. (7) and (8). We note that for 
Vectra/Thermocarb 23 composites (up to 65 vol %), C 
= 0.4638 W nT 1 K" 1 and D = 4.9256, for Vectra/ 
Fortafil 243 carbon fiber 24 composites (up to 55 vol 
%), C = 0.4841 W m" 1 K" 1 and D = 2.1478, and for 
Vectra/Panex 30 carbon fiber 24 composites (up to 55 
vol %), C = 0.4927 W nT 1 K" 1 and D = 2.7933. It is 
noted that in eq. (8) all thermal conductivities must 
have units of W m -1 K" 1 . Based upon these results, 
a new model can be developed for the in-plane ther- 
mal conductivity: 

fen = (9) 
^through 

with ^through given by K in eq. (7). 

The results of this model are shown in Figures 3 
and 4 and will now be discussed. There is a very 
good agreement between the model and the experi- 
mental data. It is noted that the models for the resins 
containing carbon fiber tend to level off as the vol- 
ume fraction is increased, indicating a possible limit- 
ing value of kfo. Mathematically, this occurs because 
the rate of growth of fehxough as predicted by eq. (7) 
exceeds that of the product ^through fen in eq. (8). This 
situation may be alleviated by choosing an alternate 
form of eq. (7c) for the through-plane thermal con- 
ductivity. 

Equations (7) and (8) are empirical models based 
upon experimental results. However, there is cur- 
rently limited through-plane experimental data avail- 
able for polymer composite systems with high load- 
ing levels, and for in-plane experimental data at any 
loading level. A rigorous parametric study investi- 
gating various systems with differing filler orienta- 
tion, length, aspect ratio, and thermal properties is 
required to obtain a better fundamental understand- 
ing of composite thermal conductivity and allow for 
development of more theoretically based models. 



We note that for carbon black the maximum vol- 
ume fraction tested was low enough that the value 
of ..<|> m = 0,64 used in previous work 31 ' 32 was not 
exceeded. However, the thermal conductivity of car- 



CONCLUSIONS 

In this project, four different carbon fillers (Kerjen- 
black carbon black, Thermocarb synthetic graphite, 
Fortafil 243 carbon fiber, and Panex 30 carbon fiber) 
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were tested for thermal conductivity at single filler 
loadings of up to 75 wt % (65 vol %) in Vectra 
A950RX LCP. Through-plane and in-plane thermal 
conductivities were measured. When comparing the 
results of resins containing Thermocarb, Fortafil 243, 
or Panex 30, the Therrhocarb composites exhibited 
higher through-plane and in-plane thermal conduc- 
tivities over the entire range of fillers studied. 

The Panex 30 composites also showed improved 
through-plane thermal conductivity when compared 
with Fortafil 243 composites at filler loadings above 
40 vol %. At lower loadings, the through-plane ther- 
mal conductivities were identical. However, the 
Panex 30 composites always had higher in-plane 
thermal conductivity. These differences may be 
attributed to the slightly higher thermal conductivity 
of Panex 30 when compared with Fortafil 243 (Tables 
IV and V). 

Two model results were presented here: the Niel- 
sen model for the through-plane thermal conductiv- 
ity and a new model for the in-plane thermal con- 
ductivity. This new model utilizes a previously 
developed correlation for the square root of the 
product of the i n-plane an d through-plane thermal 
conductivities ^/Ar^through and the Nielsen model. 
The results showed good agreement with the experi- 
mental data. In future research, these models will 
form the basis for the development of new thermal 
conductivity models for composites with multiple 
fillers. 
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also acknowledge the following undergraduate students 
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Amanda Zalud. 
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